Graphene-based heterojunction between two topological insulators 
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Quantum Hall (QH) and quantum spin Hall (QSH) phases have very different edge states and, 
when going from one phase to the other, the direction of one edge state must be reversed. We 
study this phenomenon in graphene in presence of a strong perpendicular magnetic field on top of a 
spin-orbit (SO) induced QSH phase. We show that, below the SO gap, the QSH phase is virtually 
unaffected by the presence of the magnetic field. Above the SO gap, the QH phase is restored. 
An electrostatic gate placed on top of the system allows to create a QSH-QH junction which is 
characterized by the existence of a spin-polarized chiral state, propagating along the topological 
interface. We find that such a setup naturally provides an extremely sensitive spin-polarized current 
switch. 

PACS numbers: 72.80.Vp, 72.25.Ba, 73.43.Nq, 73.22.Pr 



I. INTRODUCTION 

Electronic properties of graphene and topological in- 
sulators have received considerable attention these last 
few years [U [2]. Graphene is a two-dimensional crys- 
tal whose electronic band structure is that of a gap- 
less semiconductor, with conduction and valence bands 
touching each other at two inequivalent points K and K' 
commonly referred to as valleys. The energy of charge 
carriers vanishes at these points and disperses linearly 
with momentum in their vicinity, forming a so-called 
Dirac cone. Low-energy excitations are massless Dirac 
fermions. These exotic quasi-particles carry a topological 
Berry phase which has been shown to give rise to many 
unusual transport phenomena such as the suppression of 
backscattering (also known as Klein tunnelling) O |4], 
weak ant i- localization O |6], and a "relativistic" quan- 
tum Hall effect [ZHS]. Some of these properties, however, 
are not robust to the presence of disorder - as soon as the 
latter is sufficiently short-ranged to induce valley-mixing 
- due to the existence in graphene of an even number of 
Dirac cones. 

On the other hand, topological insulators are a newly 
discovered type of system which are insulating in the bulk 
and characterized by the existence of robust gapless ex- 
citations at their surface [2) [TOHMj. In two dimensions, 
the topologically insulating phase possesses gapless states 
propagating along its edges. There are two famous ex- 
amples of this phase: the quantum Hall (QH) insulator, 
which can be obtained by applying a strong magnetic 
field perpendicular to the plane and is characterized by 
chiral spin-degenerate edge states (Fig. lib), and the time- 
reversal-symmetric quantum spin Hall(QSH) insulator, 
which is induced by a strong spin-orbit (SO) interaction 
fiO\ \T2\ [15 and is characterized by states with opposite 
spins propagating in opposite directions (Fig. [l^). 

Whereas examples of three-dimensional topological in- 
sulators have been found to abound in Nature, two- 
dimensional systems exhibiting the QSH phase are so far 
limited to HgTe/CdTe heterostructures which only a few 
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Figure 1: Cartoon for the edge states in a 4-terminal Z-shape 
geometry when the system is in the a) QSH phase and the b) 
QH phase. In the former case, opposite spins (thick blue and 
red dashed lines) propagate in opposite directions on a given 
edge, while in the latter case they propagate in the same di- 
rection, c) Topological heterojunction, with QSH edge states 
on the left and QH edge states on the right. This junction 
can be experimentally achieved by applying a top gate (grey 
shaded region) on the right half of the sample. While one of 
the spin species (thick blue) can propagate through this junc- 
tion, the other one (red dashed) cannot and therefore gives 
rise to a chiral state, localized at the interface between QH 
and QSH phases, which connects both edges. 



experimental groups in the world can synthetize. It was 
recently suggested [16 that graphene might come to the 
rescue, much like it did with two-dimensional electron 
gases which were until a few years ago limited to epitaxi- 
ally grown semiconducting heterostructures. Graphene's 
versatile fabrication methods and insensitivity to am- 
bient temperature and chemical conditions have indeed 
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raised important expectations for all sorts of applications 
in electronics . Following the seminal observation by 
Kane and Mele that graphene with SO coupling should 
be a QSH insulator [10 , Weeks et al. numerically pre- 
dicted the possibility of opening a substantial QSH gap 
in graphene by depositing on it heavy adatoms (such as 
In or Tl) which can locally induce strong SO coupling in 
the system. In a previous paper [18 , we studied in detail 
the transport in graphene randomly covered by diluted 
adatoms and showed that the QSH phase was insensitive 
to the inhomogeneity in the coverage nad and that, for 
all purposes, the system could be mapped to a homoge- 
neous phase with a renormalized SO coupling strength 
^so = ^so^ad- We explore in this article the fate of the 
QSH phase when a strong perpendicular magnetic field 
is applied to graphene. Surprisingly, we show that the 
QSH phase is preserved for energies below the QSH gap, 
even for extremely strong magnetic fields or in the pres- 
ence of disorder. As we will see in the following, this 
can be traced back to the peculiar "relativistic" aspect 
of the QH effect in graphene (as well as to the negligible 
contribution of Zeeman splitting in graphene). The of- 
ten quoted simplified picture of SO coupling acting as an 
effective magnetic field with opposite signs for opposite 
spins clearly breaks down in this situation. Next, we turn 
to the investigation of the transport properties of a junc- 
tion between a QSH and a QH phase (Fig.jl]^). We show 
that this setup features a robust state, localized at the in- 
terface between the two topological insulators, analogous 
to the ambipolar "snake" states which arise in graphene 
quantum Hall n-p junctions [19], and take advantage of 
it by demonstrating how it can serve to realize a topo- 
logically protected spin-polarized charge-current switch. 

The paper is organized as follows. In section [Hj we 
recall the basic features of the QH and QSH phases and 
introduce the physical question at the root of this work. 



In section III, we discuss the band structure of our sys- 



tem and show in section IV how the use of topological 
Chern numbers can provide us with a clean and simple 
interpretation for it. Numerical transport calculations 
are presented in section |V| in which we show that all 
regimes (below, above, and at the topological transition) 
possess characteristic signatures. The effects of geome- 
try and disorder are discussed. Transport in the pres- 
ence of a heteroj unction between our two topologically 



insulating phases is analyzed in section [VI) We discuss 
similarities with other setups such as quantum Hall n- 
p junctions, and comment on the possibility of realizing 
a spin-polarized current switch with huge on/off ratios. 
Our results are summarized in section [VTll 



of our knowledge, however, transport signatures of the 
competition between QH and QSH phases in graphene 
have not been considered as of yet, with the exception of 
the work by Abanin et al. [24 where the QSH phase arose 
from a different mechanism (Zeeman splitting), which is 
extremely weak [42 and can therefore easily be destroyed 
by local fiuctuations of the magnetic field. In contrast, 
we study the model introduced by Kane and Mele [TO], 
to which we add the presence of a strong perpendicular 
magnetic field: 



H = VriUxCTxTz + ^yCTy) + A^oCTzTzSz 



(1) 



n = p + eA is the generalized momentum which ac- 
counts for the presence of the magnetic vector poten- 
tial A associated with a perpendicular magnetic field 
B = 5z (V X A = B), Aso is the SO-induced QSH 
gap, and vp = 3tac-c /{'^^) is the Fermi velocity (ex- 
pressed as a function of the microscopic lattice param- 
eters t (nearest-neighbor hopping amplitude) and ac-c 
(nearest-neighbor distance) which we choose in the fol- 
lowing as our working units of energy and length, re- 
spectively), {cr, r, s} are Pauli matrices in, respectively, 
sublattice, valley and spin spaces. 

The band structures of a graphene ribbon in the QSH 
and QH phases are quite different. In the QH phase, the 
perpendicular magnetic field gives ris e to La ndau levels 
Cn = ±{hvF /Ib)\/^\^\^ with Ib = \/h/ (e5), which be- 
come dispersive close to the edges of the graphene ribbon 
(Fig. [2|l). When the Fermi level is placed between two of 
these Landau levels, transport is characterized by spin- 
degenerate edge states as in Fig. [T]3, which propagate in 
a direction imposed by the sign of the magnetic field. 
In the QSH phase, the band structure consists of hyper- 
bolic bands above the QSH gap and a pair of linearly dis- 
persing ones below it (Fig. ^) . These linear bands cor- 
respond to spin-polarized states, localized at the edges 
of the graphene ribbon on a characteristic length scale 
^so = ^^f/^so' When the Fermi level is below the QSH 
gap, transport in the system can be described by counter- 
propagating edge states as shown in Fig. [l^. By compar- 
ing Figs.[l^ and[T]3, one observes that one spin species - 
hereafter referred to as the "unhappy" spin (dashed red 
in Fig. [T]) - has to reverse its direction of propagation 
when going from one phase to the other. This article 
is devoted to the study of how this reversal happens in 
two different setups: a homogeneous sample, where the 
quantum phase transition between QH and QSH phases 
is driven by electrostatic doping, and a heteroj unction 
between the two phases. 



II. HOW TO REVERSE THE DIRECTION OF 
PROPAGATION OF AN EDGE STATE ? 

The issue of how time-reversal-symmetry breaking can 
affect the QSH phase has been addressed previously in 
the literature in different settings ^2Qrf23^ . To the best 



III. BAND STRUCTURE 

The first insight is given by the spectrum of Eq. [l] 
in the presence of graphene edges. The energy spec- 
tru m for a bulk system described by Eq. IT] reads = 
±y/A%\n\ + A2^, with = {hvp /Ib)V^ T^^e lowest 
level n = stands out from the others, as each branch 
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Figure 2: Band structure of a (semi- metallic armchair) 
graphene ribbon in the a) QH and b) QSH phase. When both 
magnetic field and SO coupling are present c), the resulting 
band structure leads to a QSH phase for \Ef\ < Ago (shaded 
region) and a QH phase for \Ef\ > Ago- Compared to the 
pure QH and QSH cases, the spin degeneracy is lifted (blue 
thick and red dashed lines), which is particularly prominent 
in the lowest band which consists of spin-polarized branches 
at ^ = ±Aso- As the Fermi energy crosses the SO gap, the 
localization of the "unhappy" spin (red dashed) shifts from 
one edge to the other, while it is fully localized in the bulk 
when Ef = Ago- This is illustrated in the corresponding 
current-density plots d). On the other hand, the "happy" 
spin (thick blue) gets increasingly localized on the same edge 
as the Fermi energy crosses the transition region (not shown) . 
The parameters used are Ago = 0.02, — 8, and a ribbon 
width W = 40^. 



counter-clockwise. Note that in order for the system to 
remain gauge invariant, Peierls substitution has to be 
done on all hopping matrix elements: nearest-neighbor 
and (SO) second nearest-neighbor. To compute trans- 
port properties numerically, we make use of the software 
KNIT, which is based on an advanced recursive Green's 
function technique [25 , and works in the linear response 
regime. The numerical calculations are done with semi- 
metallic armchair boundary conditions, but our results 
are qualitatively unaffected by this choice. An impor- 
tant technical point is that the magnetic field should be 
present in the entire sample, including in the leads, in or- 
der to avoid spurious reflection at the lead-sample inter- 
face. The multi-terminal Peierls subsitution prescription 
allowing to do that is described in the Appendix. 

The full tight-binding band structure of a semi-metallic 
armchair graphene ribbon described by the Hamiltonian 
in Eq. ([2| is shown on Fig. It can be summarized 
in very simple terms: for Fermi energies inside the SO 
gap \Ef\ < Ago (shaded region), the system is in the 
QSH phase, with opposite spin channels on a given edge 
propagating in opposite directions, while for energies 
\Ef\ > Ago, the system is in the QH phase, with opposite 
spin channels on a given edge propagating in the same 
direction. Hence for a given value of Ago, the transition 
between the two phases is governed solely by the Fermi 
energy and does not depend at all on the value of the 
magnetic field (once again neglecting Zeeman splitting, 
which is very small in graphene). This quite remarkable 
result is a direct consequence of the existence in graphene 
of a ^-independent zero-energy Landau level: as soon as 
Ago 7^ 0, the spin degeneracy of the zero-energy Landau 
level is lifted, as opposed to all other Landau levels which 
remain spin degenerate [22 . This lifting leads to a QSH 
phase in the corresponding SO gap, as can be understood 
with the help of topological invariants. 



(zb) can be shown to host only one of the two spin species 
[22^ (Fig. Whereas n > levels will all disperse in 
the same direction when confinement is taken into ac- 
count, the fate of the lowest level is more subtle. To be 
more quantitative, we make use of a tight-binding model 
on the graphene hexagonal lattice, which in the presence 
of both Kane-Mele SO coupling and perpendicular mag- 
netic field can be written as 

H = -t ^ e'^^^'4,c.Cj,c.+^Ago Yl ^^je'^"4,c.^a/3C,,/3 

(2) 

Indices (z, j) label lattice sites, (a, /3) label spin indices, 
while symbols ( ) and (( )) respectively refer to nearest- 
neighbor coupling (with hopping amplitude t) and next- 
nearest-neighbor coupling (with SO-induced hopping am- 
plitude Ago = Ago/(3A/3) [10 ). The Peierls phase = 
{e/h) J^' A • dr takes into account the contribution from 
the magnetic flux threading the lattice, and = ±1 
depending on whether sites are coupled clockwise or 



IV. TOPOLOGICAL ORDER 

In this section, we discuss the topological order of the 
phases obtained by varying the chemical potential in the 
energy spectrum of Fig. In particular, we relate the 
unique transition between the QSH topological order and 
the QH topological order in presence of SO coupling to 
the specificities of the QH physics of Dirac fermions in 
graphene. Let us start by recalling the standard topo- 
logical number characterization of Landau levels when 
Ago = 0. Each Landau level n and its associated eigen- 
functions over the first Brillouin zone are characterized 
by a topological invariant, the so-called Chern number 
[26 . This topological number takes a value Q'^} = +1 
for each Landau level, independently of the Landau n, 
valley r or spin s indices. For each value of the Fermi 
energy, we can characterize the corresponding phase by a 

topological number C = ^ Zle^<£;F obtained by 
summing the Chern numbers of all filled energy bands 
[26j . For graphene and any Dirac system, however, this 
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procedure would yield an ill-defined topological number 
C due to the presence of an infinite number of filled Lan- 
dau levels below Ep. As shown recently through the use 
of non- commutative Berry's connection [27 , the correct 
topological number for a single Dirac cone takes a value 
Gr,s = -1/2 for 6-1 < Ef < eo, and Gr,s = +1/2 for 
eo < Ep < ei, where are the usual Landau levels and 

^r,s = "^er^KEp ' discussiou bclow, wc will 

make use of the more convenient topological Chern num- 
ber per spin species, Cg = Qr,s which takes values 
two-fold larger. 

Let us now turn to the energy spectrum of Fig. ^ 
where Ago 7^ 0. The presence of the SO coupling does 
not modify any of the Chern numbers per Landau level, 
but it lifts the spin degeneracy of the n = Landau 
level into the two levels dit E = iAgo- As the z compo- 
nent of spin is conserved, the topological Chern numbers 
per spin species introduced above turn out to be useful 
quantities to characterize the topological order in this 
new spectrum. They read 



-1,6; 

+i,e; 

+1,6; 
+3,6; 



-1 for e_i < Ef < -A^o , (3a) 

-1 for - Aso <Ef< Aso , (3b) 

+1 for Aso < Ef <ei , (3c) 

+3 for ei<EF <e2 , (3d) 



this time with the modified Landau levels introduced 
at the beginning of section |III| The difference between 
and C; for \Ef\ < Ago signals the appearance of a 
Z2 topological order characteristic of the QSH phase. In- 
deed, when the z component of spin is conserved, the Z2 
topological index characterizing the QSH phase is defined 
diSiy = (e^ - e;)/2 (mod 2) [2|128^. As ah Landau levels 
n ^ are still spin-degenerate, we have = C; and thus 
u = for all Fermi energies \Ef\ > Ago- For these values 
of Ef^ the system lies in a QH phase characterized by the 
usual topological Chern number C = + C;. However 
for \Ef\ < Ago, Eq. Q leads to a non- trivial Z2 in- 
dex = 1, while the total Chern number simultaneously 
vanishes C = + C; = 0. The system then lies in a dif- 
ferent topologically insulating phase: the QSH insulator. 
This shows that as the Fermi energy crosses the values 
±Aso, the system undergoes a quantum phase transition 
between two topological insulators: a QH phase and a 
QSH phase: 

u = 0,e = -2 for e_i < Ef < -A^o : QH , (4a) 
z/ = l,e = for -Aso<^F<Aso : QSH , (4b) 
= 0, e = 2 for Aso <Ef <ei : QH , (4c) 
u = 0,e = 6 for ci < Ef < €2 : QH . (4d) 

This transition appears crucially tied to the Dirac physics 
of graphene and the presence of the n = Landau level: 
in the present case, we do not need a SO coupling to 
overcome an energy gap in order to drive this transition, 
as would be the case for non-relativistic fermions with 
Landau gap or if graphene had a trivial mass gap 
maz. The spin degeneracy lifting of the n = level is all 



that is required here. Let us note finally that while this 
argument formally uses the conservation of the Sz spin 
component, the robustness of topological numbers proves 
it to remain valid if non-S'^-conserving terms are included 
in the Hamiltonian, provided two branches of the n = 
Landau level with opposite spins remain non-degenerate. 



V. TRANSPORT SIGNATURES 
A. Ballistic regime 

Let us now study how this topological phase transition 
appears in transport. We will study the multi-terminal 
(dimensionless) differential conductance Tab which ex- 
presses how much current dia is collected in lead a when 
the voltage in lead b is raised by dVb, 



(5) 



Additionally, in order to observe the edge states directly, 
we also study (in color plots) the differential local cur- 
rent density di{f)/dVa which allows to clearly observe the 
edge states inside the sample. In Fig. |2]i, we show local 
current-density plots which illustrate how the behavior 
of the "unhappy" spin changes as a function of the Fermi 
energy: it goes from propagating along one edge when 
Ef < Ago (left panel) to propagating along the other 
as Ef > Ago (right panel), while it gets localized in the 
bulk at the critical point Ef = Ago (middle panel). We 
find that the energy window where this localization is 
observed is extremely narrow and decays exponentially 
with the width of the sample. Note that this scenario is 
completely different from what one would expect start- 
ing from the naive toy model of SO coupling acting as a 
spin-dependent magnetic field B^oZSz. In this case, for 
a critical value of the real magnetic field B = B^o^ the 
"unhappy" spin would feel no magnetic field at all and 
be fully delocalized. On the contrary, we observe that 
the QSH phase is virtually independent from the real 
magnetic field and that the "unhappy" spin actually gets 
localized when Ef = Ago, illustrating the limitations of 
the toy model in this situation. 

In the vicinity of the transition, the "unhappy" spin 
keeps propagating along a given edge but its classical 
cyclotron orbit center Xc = —kyl^ (with ky the longitu- 
dinal wave vector component) is shifted inwards as the 
Fermi energy increases (see Fig.[2]3). On the other hand, 
the "happy" spin gets increasingly localized on the same 
edge when the Fermi energy increases, which is qualita- 
tively equivalent to the usual QH case. One can indeed 
show that the notion of a classical cyclotron orbit cen- 
ter remains well defined here, despite the presence of SO 
coupling, and that the corresponding eigenstates are very 
similar to those found in the pure QH regime [22 • 

Before presenting the rest of our numerical data on 
transport in the vicinity of the QSH-QH transition, a 
few words should be said about the shape of the spectrum 
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corresponding to the non-zero Landau levels. While the 
branches corresponding to the "happy" spin seem basi- 
cally unaffected by the SO coupling, the two branches 
of the "unhappy" spin display very different behaviors 
(Fig. [2]3). In particular, one of them is significantly bent 
by the SO coupling, such that counter-propagating states 
along the same edge appear in a finite window of energy. 
This leads to the possibility of backscattering and there- 
fore destroys the robustness of the QH phase in this en- 
ergy window, the size of which can nevertheless be signif- 
icantly reduced by increasing the width of the graphene 
ribbon. Note that at negative energies, Ef < —Ago, 
"happy" and "unhappy" spin species exchange their roles 
(see Fig. (2)3). 

We proceed to investigate further how the transition 
between QSH and QH phases shows up in transport, 
and consider the 4-terminal ^-shape geometry depicted 
in the inset of Fig. [3^. In the core of Fig. |3^, we plot 
as a function of the Fermi energy the current collected 
in leads 1, 2 and 3 when injected from lead 0. Noth- 
ing unexpected happens away from the transition: the 
transmission coefficients feature characteristic signatures 
of current-carrying edge states (left and right panels of 
Fig. [2]i). Around the critical value Ep = Ago, however, 
we observe in Fig. [3^ that the "unhappy" spin is fully 
transmitted in lead 2. This can be understood as fol- 
lows. As the transition point is approached, the classical 
cyclotron orbit center is shifted (inwards) away from the 
edge (as discussed above) by a distance which can reach 
2Rc, where Rc = l%/£,so is the classical cyclotron radius 
of the n = level. When 2Rc > Wi/2, where Wi is 
the width of lead 1, then, somewhere in the vicinity of 
the transition, the incoming state cannot penetrate lead 
1, hence leading to perfect direct transmission to lead 2 
as observed in Fig. [3^. This feature could in principle be 
checked experimentally by simply varying the strength of 
the magnetic field close to the transition: for wider sys- 
tems or at larger magnetic field, the width of the peak 
of direct transmission from to 2 shrinks and eventually 
disappears. 

A clear way to distinguish between QH and QSH 
phases in a transport measurement is to plot the "di- 
rect" dimensionless differential conductance from lead 
to lead 2 as a function of the classical cyclotron radius 
Rc oc B-^. This is illustrated in Fig.^jo: below the QSH 
gap, the "direct" transmission is zero for any value of the 
magnetic field; on the other hand, above the QSH gap, 
transmission can become non-zero as the magnetic field 
weakens, due to the breakdown of the QH effect when Rc 
becomes larger than the width of the system. It thus in- 
terpolates between zero (in the QH regime) and a finite 
value which depends on the number of occupied bands 
at the Fermi level in the ribbon when B = 0. Of course, 
in the latter limit, the topological protection is lost and 
the value of the transmission will strongly depend on the 
disorder configuration. 




Figure 3: QSH-QH transition in a 4 terminal ^-shape ge- 
ometry, a) Dimensionless differential conductance from lead 
0, where the current is injected, to outgoing leads 1, 2 and 
3 in the geometry depicted in the inset, as a function of 
the Fermi energy (Ago = 0.02, Ib — 8, and lead widths 
Wo = 64^3, Wi = 22\/3). Below the SO gap, current is car- 
ried by QSH edge states, which propagate on different edges 
(left inset), while above the SO gap it is carried by QH edge 
states, which propagate on the same edge (right inset), b) 
Dimensionless differential conductance from lead to lead 2 
as a function of the classical cyclotron radius Rc (same val- 
ues of Ago, Wo and Wi as in a). Circular symbols (thick 
line) correspond to data with Ep < Ago, while other symbols 
(dashed lines) correspond to data with Ep > Ago, respec- 
tively Ef/^so — 1-06 (blue diamonds), Ep/^so — 1-25 (vio- 
let triangles), and Ep/^so = 1-44 (red squares). While this 
"direct" transmission is vanishingly small and independent of 
the magnetic field below the SO gap, it increases with both 
Rc and Ep above the SO gap. The latter situation arises be- 
cause the QH phase is destroyed as 2Rc > VKi, and because 
the number of transmitting channels increases with Ib and 
Ef- 



B. Effect of disorder 

In this section, we study the robustness of the results 
described in the ballistic regime with respect to the pres- 
ence of various types of disorder. In the absence of a 
magnetic field, the QSH phase is known to be very ro- 
bust with respect to the presence of scalar disorder (i.e. 
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disorder that breaks neither time reversal symmetry nor 
spin rotational symmetry), as introducing backscatter- 
ing between edge states involves tunneling through the 
gapped bulk region. On the other hand, disorder that 
breaks both time reversal symmetry and spin rotational 
symmetry (such as magnetic impurities) leads to scatter- 
ing between the two counter-propagating edge states of 
the same edge, which leads to intra-edge backscattering. 
If this intra-edge backscattering becomes strong (or the 
edges are very long), edge states may eventually get lo- 
calized, such that no edge transport occurs anymore. In 
this paper, as we explicitly break time reversal symmetry 
with a magnetic field, disorder that breaks only spin ro- 
tational symmetry (such as Rashba-like SO terms arising 
from adatoms) effectively behaves as magnetic impurities 
and could potentially lead to the same breakdown of the 
QSH phase. 

In order to study these effects quantitatively, we add 
the general perturbation 

ifdis=E4j E (6) 

i,Q!,/3 \fi=0,x,y,z / 

to our tight-binding Hamiltonian ([2|. The onsite po- 
tentials Vi^y, are independent random variables uni- 
formly distributed inside a given interval on each 
site of the system (Anderson disorder). We study 
three different sorts of disorder with different symme- 
tries: l^^o ^ [— 1^/2] takes into account scalar 
(spin-independent) disorder, Vi^^ G [— V2/2,14/2] rep- 
resents Zeeman-like (Sz conserving) impurities, and 
{Yi.xi ^i,y\ ^ [— Kn/2, Kn/2] capturcs the infiuence of 
non-conserving impurities. 

The results are presented in Fig. |4] for the 4-terminal 
^-shape geometry (depicted in the inset of Fig. [3^). Note 
that in order to obtain a significant effect of disorder, we 
used extremely (unrealistic) high values of disorder, much 
higher than the SO gap itself. A first general qualitative 
conclusion is that the ballistic results presented above 
are extremely robust with respect to all kinds of disor- 
der. More precisely, we find, as expected, that extremely 
large values of Vs or are needed to significantly affect 
the transport coefficients. "Magnetic" impurities {Vm)^ 
which can induce intra-edge backscattering, do affect the 
QSH phase at weaker values. We note, however, that 
the system remains topological (i.e. does not become a 
simple ordinary insulator) as transport is still dominated 
by edge contributions: for instance, at Vm = 0.5, the 
direct transmission probability T20 (from lead to lead 
2) remains many orders of magnitude smaller than Tio 
(from lead to 1). The topological nature of the phase 
is encoded in the Z2 topological invariant, which cannot 
change unless the bulk gap is closed by the perturbations 
we consider. We checked (not shown) that the above re- 
sults are essentially unaffected upon changing the SO gap 
by a factor of 2, putting all three sorts of disorder simul- 
taneously or changing the Fermi energy (inside the SO 
gap). 




Figure 4: Effect of very strong disorder on the QSH phase in 
the 4-terminal ^-shape geometry. We plot the dimensionless 
differential conductance from lead to leads 1 and 2 as a func- 
tion of disorder strength for Ago = 0.02, /b — 8, Ep = 0.05, 
and otherwise same system as in Fig. [3] Circular symbols 
stand for Tio and triangular symbols for T20. The behavior for 
T30 (not shown) is the same as for Tio- Dashed (black) lines 
correspond to purely scalar disorder (Vs), dotted (red) lines 
correspond to purely Zeeman-like disorder (Vz), and filled 
(blue) symbols correspond to purely Sz non-conserving dis- 
order (Vm)- The latter is clearly the dominant effect on Tio, 
though unrealistically large values are required to quantita- 
tively affect the edge state transmission (see main text for a 
mean free path estimate) . Each data point has been averaged 
over 48 disorder configurations. Unless shown, error bars are 
smaller than symbol sizes. Full (black) line: fit ^ = 1 — S.Ox^. 



To rule out any concerns raised by the intra-edge scat- 
tering due to Vm^ the corresponding intra-edge mean 
free path (or equivalently the localization length, as 
both are roughly equal for one-dimensional states) can 
be estimated using Fermi's golden rule: /g = ^dTe 
where Vd = h~^{dE/dk) is the drift velocity of the 
edge states (extracted for instance from Fig. [2|, and 
h/re (X {Vm)'^\dk/dE\. Alternatively, one may extract it 
directly from the numerical calculations of Fig. |4j since 
Tio ^ ^ — {L/le) (with L the length along the edge, 
roughly 5 nm in this instance). Fitting the small Vm 
regime with Tio ^ 1 — 3.01^, we obtain (in nm) 

le « (KnA)"' , (7) 

where we have explicitly restored the hopping amplitude 
t (our energy unit) in order to get numbers. The ballis- 
tic results are essentially unaffected for systems smaller 
than /g, while the edge states become localized for larger 
systems. In the context of a QSH phase induced by 
adatoms (using for instance the indium atoms proposed 
by Weeks et al. [E]), a possible source of Frn-hke disor- 
der comes from the SO coupling induced by the adatoms 
themselves. Typical values for Vm are smaller than 1 
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meV, which translates into extremely large intra-edge 
mean free paths 1^ > 1 mm. Hence we estimate that 
this perturbation should be largely irrelevant in real size 
samples. 

To summarize, all of our results are essentially unaf- 
fected by the presence of disorder, except when the Fermi 
energy lies in the vicinity of the QSH-QH transition, in 
which case strong disorder can give rise to a random net- 
work of QH and QSH regions through which a percolat- 
ing cluster connecting opposite edges can therefore lead 
to backscattering [43^ . 



VI. TOPOLOGICAL HETEROJUNCTION 

We take advantage of the above described topologi- 
cal quantum phase transition as a function of the Fermi 
energy to propose a setup which allows for a direct junc- 
tion between two different topological phases in the same 
sample. Let us consider the case in which an additional 
electrostatic gate enables to split the system in two parts, 
one in which the Fermi level is in the QSH phase, and 
the other in which the Fermi level is in the QH phase 
(Fig. [l]^). This constitutes a QSH-QH junction which 
shares some similarities with quantum Hall n-p junctions 
previously fabricated in graphene [29l |30]. Indeed, the 
incoming "unhappy" spin at the junction has no other 
choice but to propagate along the interface in order to 
reach the only other available channels which lie on the 
opposite edge. This is reminiscent of the situation both 
spin channels have to face in the QH regime when they 
must cross a n-p junction, since their direction of prop- 
agation on a given edge is reversed for negative ener- 
gies. Various theoretical models have been proposed in 
the latter setup [3Tti34] but they all fail to explain the 
experimental observations [29l [30] , probably due to some 
dephasing mechanism taking place in the vicinity of the 
Dirac point which is obscured by charge-density fluctua- 
tions (so-called electron- hole puddles). The system we 
consider could therefore provide a new perspective to 
solve this puzzle, as the QSH-QH transition takes place 
at a value of energy which can be far away from the Dirac 
point [44] for realistic values of SO-induced QSH gap. 

More generally, our proposal offers the possibility of 
studying the nature of the state that propagates at the 
interface between QSH and QH phases, which are charac- 
terized by different topological invariants [TTJ [26] . What 
we usually refer to as QSH (or QH) edge states are states 
propagating between QSH (or QH) insulators and a triv- 
ial insulator (the vacuum, typically). As the QSH insu- 
lator is characterized by a Z2 number, there is only one 
QSH topological phase: junctions between QSH phases 
with different Fermi energies (including n-p junctions) 
have no effect on transport, as the spin-polarized states 
can propagate through these junctions. On the other 
hand, the QH topological invariant is a Z number, which 
counts the number of edge channels, and the notion of 
QH junctions therefore makes sense. In this case, one 




0.5 1 1.5 2 

Vg/(A3o-Ep) 

Figure 5: Heterostructure in a 4 terminal Z-shape sample 
as depicted in Fig. [l]:. Transmission probabilities from lead 
0, where the current is injected, to outgoing leads 1 (dotted 
lines), 2 (dashed lines) and 3 (full lines) as a function of the 
top gate voltage Vg. When Vg < (Ago — Ep), such that left 
and right regions are in the QSH phase, current is perfectly 
transmitted by the QSH edge states, as shown in the current- 
density plot in the left inset. When Vg is high enough (Vg > 
(Ago — Ef)) that the right part of the sample enters the QH 
phase, a QSH-QH junction is created, which is characterized 
by a chiral state propagating along the interface. This is 
illustrated in the current-density plot shown in the right inset. 
Once it has reached the opposite edge, this chiral state is 
partially transmitted in lead 2, partially transmitted in lead 3, 
with proportions shown in the main plot. The light red curves 
correspond to an abrupt voltage change across the junction 
region while the dark blue curves correspond to a smooth 
transition. Parameters are Ep/Aso = 0.58 (in the left half of 
the sample), Aso 0.02, ^ 8, and widths Wo = Wc = 40^/3 
for the leads and central region. 



expects the existence of chiral propagating states, local- 
ized at the interface corresponding to the Landau level 
crossing. For QH n-n' junctions, these states are "bub- 
bling" states [34], which simply follow the drifting Hall 
motion of charge carriers subjected to crossed electric 
and magnetic fields. For QH n-p junctions, these states 
are ambipolar "snake" states [19], which can be seen as 
classical skipping orbits of mixed electron- "hole" charac- 
ter. The characteristics of the state propagating at the 
interface between QSH and QH insulators, on the other 
hand, are still unclear as far as we know. 

Our proposal is the following. Consider the Z-shape 
geometry shown in Fig. [ij^. It is split into two regions: 
the left half and the right half. A top gate is applied 
on the latter in order to tune the Fermi level in both re- 
gions independently. In Fig. [5| we plot, for a given value 
of the Fermi energy in the left half inside the QSH gap, 
the differential conductance T^o from lead to leads 1, 
2 and 3 as a function of the gate voltage (which deter- 
mines the value of the Fermi energy in the right half). 



8 



While the Fermi energy in the right half remains below 
the value of Ago, transport characteristics remain unaf- 
fected by the gate (left inset of Fig. [5|. However, as 
soon as the Fermi energy in the right half crosses the 
QSH gap, it gives rise to a QSH-QH junction charac- 
terized by a topological state at its interface (see right 
inset of Fig. [5|. This chiral state propagates along the 
interface until it reaches the opposite edge, and then gets 
partially transmitted in lead 2, partially transmitted in 
lead 3, with proportions which depend on the micro- 
scopic details of the model (Fermi energies, the length 
of the interface, the smoothness of the potential step, 
the amount of disorder, the possible valley- space polar- 
ization of the edge states, etc), the study of which is 
left to subsequent work. This system constitutes a very 
efficient spin-polarized charge-current switching mecha- 
nism, as the current in lead 1 (respectively 2) can be 
reversibly switched from one (respectively zero) to zero 
(respectively non-zero) while simultaneously being spin- 
polarized (see Fig. |5|. Additionally, this switching can 
be activated by simply tuning the voltage in the top gate 
over a very small energy range (whose lower bound will 
depend on the magnitude of the disorder in the vicinity of 
the transition), and should yield extremely sharp transi- 
tions with values of on/off current topologically protected 
from the presence of disorder |45j . 



VII. CONCLUSION 

In summary, we showed that the transition between 
QSH and QH phases in graphene is independent of the 
value of the magnetic field (neglecting the weak effect of 
Zeeman splitting) and can be crossed simply by tuning 
the value of the Fermi energy across the SO-induced QSH 
gap. In a heteroj unction, one of the spin species gives rise 
to a chiral state propagating along the interface between 
QSH and QH phases. The nature of this special state is a 
fascinating issue which could bring new light concerning 
the unresolved mystery of conductance plateaus observed 
in quantum Hall n-p junctions [29^, ^30^ . We also showed 
that the tunable transition between the QSH and QH 
topological phases could serve as a spin-polarized charge- 
current switch with potentially extremely high, topolog- 
ically protected, on/off ratios. An interesting future di- 
rection of research could be to investigate whether this 
tunable topological phase transition can arise in bilayer 
graphene [35], which also possesses zero-energy Landau 
levels [36^ and has very recently been shown to host a 
Z2 topologically insulating phase ^ for strong enough 
Rashba SO coupling. 



Appendix: magnetic field in multi-terminal 
calculations 



In this appendix, we explain how to compute the 
Peierls phase (pij between two atoms i and j in multi- 



terminal systems [38]. It is given by 

A • dr 



(8) 



where = (xi^i/i) is the spatial position of site i. A 
common choice for A is the Landau gauge (Lg) 



ALg{x,y) = -By 



which leads to 



-B{xi — Xj) 



Vi + Vj 



(9) 



(10) 



(using linear paths between atoms). The numerical pre- 
scription follows simply: one calculates the real coordi- 
nates of the two atoms and uses the above equation to 
get . An important aspect of the Landau gauge is that 
the phase depends on the x-coordinate only through the 
difference of x between the two atoms i and j. This 
is crucial for taking magnetic field into account in the 
leads: the leads are semi-infinite periodic systems made 
of layers. They are described by an intra-layer Hamilto- 
nian Hq and an inter-layer Hamiltonian V . Within the 
Landau gauge, we find that the matrices Hq and V of 
horizontal leads are layer-independent, which allows the 
use of standard schemes to calculate their self-energies. 
However, general samples (such as the ^-shaped sample 
studied in this article) can have leads with an arbitrary 
angle with respect to the ?/-axis. For those leads, the 
corresponding gauge choice is 



A(r) = -5(r.e2)ei 



with 



ei 



cost 
sin^ 



62 



— smt 
cos6> 



which leads to 



^3 



where the potential 

<I>i = 5(l-cos2l9) 



B sin 26*- 



(11) 



(12) 



(13) 



(14) 



is a pure gauge potential allowing to go from one choice 
of gauge to the other. A (possible) general prescription 
for an arbitrary system is now the following: one assigns 
a potential <I>^ = to all sites except those belonging 
to a lead. For the latter, one uses Eq. ( [l4| ) with the 
appropriate angle 0. Then one calculates the phases 
according to Eq. (13). 



All the prescriptions above are given in real space. It 
is of course possible to calculate analytically the equiva- 
lent prescriptions in terms of the integer coordinates on 
the Br avals sublattices, as both are in one to one corre- 
spondance. However, for numerical purposes, it is more 
convenient to calculate the real spac e position s of the 
atoms numerically and then use Eqs. ( 1Q|14|13" ) in order 
to use a lattice-independent prescription. 
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